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The Renyi statistics in the canonical and microcanonical ensembles is examined in the general 
case and in particular for the ideal gas. In the microcanonical ensemble the Renyi statistics is 
equivalent with the Boltzmann-Gibbs statistics. By the exact analytical results for the ideal gas, 
it is shown that in the canonical ensemble in the thermodynamic limit the Renyi statistics is also 
equivalent with the Boltzmann-Gibbs statistics. Furthermore it satisfies the requirements of the 
equilibrium thermodynamics, i.e. the thermodynamical potential of the statistical ensemble is a 
homogeneous function of degree 1 of its extensive variables of state. We conclude that the Renyi 
statistics duplicates the thermodynamical relations stemming from the Boltzmann-Gibbs statistics 
in the thermodynamical limit. 
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I. INTRODUCTION 

The general aim of this paper is to establish the 
connection of the Renyi statistics in the canonical and 
microcanonical ensembles with the equilibrium ther- 
modynamics and to compare it with the Boltzmann- 
Gibbs statistical mechanics. The Renyi statistics is 
defined by the Renyi statistical entropy [H, H[ with 
| respect to the usual norm equation for probabilities 
and the usual expectation values for dynamical vari- 
ables [H, 0. It differs from the Boltzmann-Gibbs 
statistics only in the definition of the statistical en- 
tropy. The main advantage of the Renyi statistics is 
the power-law distribution function instead of the ex- 
ponential distribution function of the Gibbs statistics. 
Therefore the Renyi entropy [3j, |J] has been applied in 
several situations, for example, fractal and multifrac- 
tal systems @, 0, high energy particle production Q , 
nuclear physics [9( and black holes [10|. However in 
these fields the Tsallis statistics is more significant due 
to the nonextensive nature of the Tsallis statistical en- 
tropy 00. 

The Renyi statistics as the theory of the equilibrium 
statistical mechanics is defined at the equilibrium in 
the maximum entropy point. Therefore it must satisfy 
all requirements of the equilibrium thermodynamics 
and must be in agreement with the theory of prob- 
ability, and some general physical principles. The 
Renyi statistics is thermodynamically self-consistent 
whenever it satisfies in the thermodynamic limit the 
zero, first, second, and third laws of thermodynam- 
ics, the principle of additivity, the fundamental equa- 
tion of thermodynamics, the Gibbs-Duhem relation, 
the differential equation for thermodynamical poten- 
tial, the Legendre transform and the Euler theorem. 
It can be proved that the zeroth law of thermody- 
namics, the principle of additivity, the Euler theorem, 
and the Gibbs-Duhem relation are valid if the ther- 
modynamical potential of the statistical ensemble is 
a homogeneous function of first degree of the exten- 
sive variables of state [Til . [2|| . For finite systems the 
thermodynamic potential is an inhomogeneous func- 
tion. Therefore finite systems are not thermodynam- 



ically self-consistent. In the thermodynamical limit 
the thermodynamical potential of the statistical en- 
semble is a homogeneous function of the first degree. 
Therefore for infinite systems the connection of the 
statistical mechanics with the equilibrium thermody- 
namics is provided. The equilibrium thermodynam- 
ics is a macroscopic phenomenological theory which is 
defined in the thermodynamic limit, while boundary 
effects can be neglected [131 ]. 

The Renyi statistical entropy is an extensive func- 
tion if it is a multivariate function of the set of proba- 
bilities. Therefore the Renyi statistics was considered 
thermodynamically self-consistent without any fur- 
ther verification. Otherwise the Tsallis statistics was 
thoroughly investigated. The studies done by many 
authors fail to answer the question whether the zeroth 
law of thermodynamics and the principle of additivity 
for the Tsallis statistics are valid. These unsuccessful 
attempts are based on the following physical concepts: 
the physical temperature [1 7]| . the entropic parameter 
q as universal constant [la, |l6| , the temperature de- 
pendent nonextensive Hamiltonian 18], the transfor- 
mations from the Tsallis entropy to the Renyi (Gibbs) 
one @, [lj| l20l. l2l| and the inconsequent thermody- 
namic limit [22l| . In 0, [2||, the author proved that 
the zeroth law of thermodynamics in the thermody- 
namic limit for the Tsallis statistics is satisfied only if 
the parameter z = l/(q — 1) is an extensive variable 
of state. Notice that the consistent thermodynamic 
limit for the Tsallis statistics was firstly proposed in 
the particular case by Botet et al. in [2J, UM ■ 

In our previous paper Q , we obtained the main re- 
lations for the Renyi statistics in the canonical ensem- 
ble in the general case and for the ideal gas. It was 
shown that the Renyi statistics in the microcanoni- 
cal ensemble is equivalent with the usual Boltzmann- 
Gibbs. In the present paper we will show that the 
Renyi statistics in the canonical and microcanonical 
ensembles in the thermodynamical limit is thermody- 
namically self-consistent whenever it coincides with 
the usual Boltzmann-Gibbs statistics. 

The authors in (2(| [StJ erroneously claimed that the 
Renyi statistics has the exponential distribution func- 
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tion instead of the power-law probability distribution 
function. In (2(| this follows as a wrong interpreta- 
tion for the distribution function of the microcanon- 
ical ensemble and in [27| the exponential stationary 
distribution function was obtained on the basis of a 
generalized entropy maximization procedure, distinct 
from the usual Renyi one. 

This paper is organized as follows. In the second 
section we briefly outline the canonical ensemble for 
the Renyi statistics in the general form and for the 
ideal gas and define the zeroth law of thermodynam- 
ics. The same results for the microcanonical ensemble 
are given in the third section of the present paper. 



II. CANONICAL ENSEMBLE (T,V, z, N) 

A. General formalism 

The equilibrium statistical mechanics is said to be 
the Renyi statistics if the equilibrium phase space dis- 
tribution function satisfies conditions imposed by the 
Renyi's statistical entropy [1, Q| with respect to the 
norm equation for a phase space distribution func- 
tion and an expectation value for a dynamical variable 
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where g is the phase space distribution function, dT = 
dxdp is an infinitesimal element of phase space, k is 
the Boltzmann constant, geRis the real parameter, 
q G [0, oo], and z = is a variable of state. 

Consider the canonical ensemble of the classical dy- 
namical systems of N particles at a constant temper- 
ature T, in a volume V, and with a thermodynamic 
coordinate z, in a thermal contact with a heat bath. 
To express the equilibrium phase space distribution 
function from the macroscopic variables of state, we 
consider the thermodynamical method based on the 
fundamental equation of thermodynamics [lil |25| in- 
stead of the Jaynes principle The fundamental 
equation of thermodynamics at the fixed values of 
T, V, z, N can be written as 



{TdS - d(H)) T ,v,z,N = 0. 
Using Eqs. fTJ and ©, we obtain 
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where the expectation value of the Hamiltonian H is 
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Since in the canonical ensemble differentials dH 
dz = dg = 0, we obtain 
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where a\,a.i and a 3 are certain constants. Substi- 
tuting Eqs. {l} and ([6]) in the first equation ([7]), we 
get 
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Averaging Eq. ([5]) and eliminating the constant at, 
we obtain the phase space distribution function 
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The norm equations pTJ)) . (fTTj) arc derived from 
Eqs. ([2]) and ([6]). Therefore the phase space distri- 
bution function is determined by two norm functions 
of the variables of state, Zr = Zr(T,V,z,N) and 
(H) = (H)(T,V,z,N), respectively. They are two 
solutions of the system of equations (fTOj) . (fTTj) . The 
thermodynamical system is totally determined by the 
free energy. It is the Legendre transform of energy 
with respect to the entropy of the system and it can 
be written as 



F = (H) - TS 
= -fcTln 

where 

S = fcln 
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The expectation value of a dynamical variable A can 
be written as 
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The fundamental equation of thermodynamics for 
Renyi statistics is derived at constant values of the 
variables of state (T, V, z, N) following the procedure 
given in Ref. [25|, we obtain 



TdS = d(H) + pdV + Xdz - fidN, 



where 
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Here the property of the Hamilton function, 
(dH/dT) v . z . N = (dH/dz) T ,v,N = 0, is used. 

The thermodynamic potential F of the canonical 
ensemble (T, V,z,N), i.e., the Helmholtz free energy 
(fT2|) . satisfies the differential relation which can be 
written as 



dF = -SdT -pdV - Xdz + fjdN, (19) 



where 



x = "O t „- *=(§) TO <2i) 

The free energy is the thermodynamic potential and 
it determines the ensemble averages of the canonical 
ensemble. 

The fundamental equation of thermodynamics (| 15(1 
provides the first and second laws of thermodynamics 

SQ = TdS, SQ = d(H) + pdV + Xdz - fidN, (22) 

where SQ is a heat transfer by the system to the envi- 
ronment during a quasistatic transition of the system 
from one equilibrium state to another nearby one. The 
heat capacity, SQ = CdT, in the canonical ensemble 
at the fixed values of T, V, z, N can be written as 
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Note that the thermodynamical relations for quantum 
Renyi statistics are the same as for the classical one 
given above. 

Let us prove the zeroth law of thermodynamics for 
the Renyi statistics in the canonical ensemble. Divide 
the system into two parts, H = Hi + H 2 , under the 
conditions 

T = Ti = T 2 , N = JVi + N 2 , V = Vx + V 2 , (24) 
z = z\ + z 2 z-extensive, (25) 
z = z\ = z 2 z-intensive. (26) 

The heat transfer SQ is additive. Therefore, from the 
second law of thermodynamics (f2"2"| it follows that the 
zeroth law, T = T\ = T 2 , is valid whenever the total 
differential of entropy dS is additive 

SQ = SQi + SQ 2 = TidSi + T 2 dS 2 = TdS, (27) 
dS = dSx + dS 2l T = T X =T 2 . (28) 

Let us show that Eq. ([28|) is valid only if the entropy 
S is an additive function 

S{T,V,z,N) = Sx(T 1 ,V 1 ,z 1 ,N 1 )+S 2 {T 2 ,V2,Z2,N 2 ). 

(29) 
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Differentiating Eq. (|29"]) with respect to variables of 
state T, V N, we obtain 



T = Ti = T 2 , (30) 
V = V x + V 2 , (31) 
N = N 1 +N 2 . (32) 



Differentiating Eq. (f2"5|) with respect to variable of 
state z, we have two possibilities 



z — z\ + z 2 , (33) 
z = z\ = z 2 . (34) 



Then from Eqs. l[3"0 ]) -(|3"g ]) and the differentials of 
Eqs. (|M|l-([26l). we get 
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Therefore the zeroth law of thermodynamics, T = 
Ti = T 2 , is satisfied whenever the entropy is additive, 
S = S 1 + S 2 . 

Let us show that the additivity of entropy (|29|) cor- 
responds to the homogeneity of the function of en- 
tropy. Let S be the homogeneous function of first 
degree of its arguments (T, V, z, N). Then for z ex- 
tensive or z intensive, we have, respectively, 

S(T, XV, Xz, XN) = XS{T,V,z,N), A = (36) 



S(T, XV, z, XN) = XS{T,V,z,N). 



(37) 



Divide the system into two parts under the conditions 
pi |) -([26 l) . From Eq. ([24]) for N, we have 



(38) 
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Then the intensive specific variables correspond to the 
additive (extensive) variables 

Ay = AiVi = X 2 V 2 , Xz = X lZl = X 2 z 2 , (39) 

where Ai = l/N x and X 2 = 1/N 2 , while A = 1/N. 
Here z is an extensive variable. Under the conditions 
of Eqs. 42§D and ((Ml)-®, it follows that the specific 
entropy s — XS is an intensive function 

XS{T,V,z,N) = X 1 S 1 (T 1 ,V 1 ,z u N 1 ) = 

= X 2 S 2 (T 2 ,V 2 ,z 2 ,N 2 ). (40) 



Conversely, from Eqs. (gDJ) and ([5S ]) -([5g ]l . we obtain 
the additivity of entropy ([29]) . Therefore the additiv- 
ity of entropy (|29|) is a consequence of its homogene- 
ity ([5tJ|) . (|3T)l (extensivity). Thus the zeroth law of 
thermodynamics for Renyi statistics in the canonical 
ensemble is valid whenever the entropy (or the ther- 
modynamic potential) is a homogeneous function of 
first degree of all extensive variables of state [3, HE] • 
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If the entropy is a homogeneous function of degree 
one (|3"6")) . (157)) , then S satisfies the Euler theorem and 
the Gibbs-Duhem relation. Differentiating Eqs. ([36]) . 
(1371) with respect to A, and putting A = 1, we obtain 
the well-known Euler theorem for the homogeneous 
functions: 



Tr 8S OS dS 
V dV +Z 0l + N 0N 
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Using the Legendre transform (Tj"2")) . the homogeneity 
property of energy (H) , and the thermodynamic rela- 
tions (pp|). (f2"Tj) . we get the Euler theorem [29[: 

TS = (H) + pV + Xz - fxN z-extensive,(43) 
TS = (H) +pV - fiN z-intensive. (44) 

Applying the differential operator with respect to the 
ensemble variables (T, V, z, N) on Eqs. (|33"]). (|31)) and 
using the fundamental equation of thermodynamics 
(IT5|) . we obtain the Gibbs-Duhem relation |29j 

SdT = Vdp + zdX - Ndfi z-extensive, (45) 
SdT = Vdp - Xdz - Ndfj, z-intensive. (46) 

Therefore in the canonical ensemble the Euler the- 
orem results from the homogeneity property of the 
entropy (pil)j) . (|3T)l . the energy (H), and the thermo- 
dynamic potential. For z being extensive the Euler 
theorem (|43[) contains the term Xz and it is con- 
sistent with the fundamental equation of thermody- 
namics ([15)) and with the Gibbs-Duhem relation (|4"5"j) . 
However for z being intensive the Gibbs-Duhem rela- 
tion (|46|) is violated, and the Euler theorem (|44f is not 
consistent with the fundamental equation of thermo- 
dynamics (Ti"5"|) . Thus the variable of state z can only 
be extensive. 

Dividing the system into two parts under the con- 
ditions (|24I26[) . i.e. assuming that z is intensive, we 
obtain that the phase space distribution function ©- 
(fTTj) does not factorize and the Renyi entropy ([I]) , (p~3|) 
is nonadditive 
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Thus, the zeroth law of thermodynamics and the prin- 
ciple of additivity for the Renyi statistics in the canon- 
ical ensemble are violated. In the general case for 
finite values of the variables of state V, z, N the ther- 
modynamical potential (fT2"]) of the canonical ensemble 
is not a homogenous function of first order. Therefore 
the Euler theorem for finite systems in the canonical 
ensemble is not satisfied, 

TS ^ (H) +pV + Xz- (j,N z-extensive, (48) 
TS ^ (H) +pV- fiN z-intensive. (49) 

The proof of inequalities P7| ~ P9")) is given below for 
the ideal gas. 

In the general case, the Renyi entropy (fT]), (p~3|) 
and the thermodynamical potential (Ti~2"|) of the canon- 
ical ensemble is not a homogeneous function. But 



it can become homogeneous in the thermodynamic 
limit. A limiting statistical procedure is called the 
thermodynamic limit if the following conditions hold: 
(i) N — > oo, V — + oo, z — * oo and v — V/N — const, 
z = z/N = const for z extensive; (ii) N — * oo, V — > oo 
and v = V/N = const, z = const for z intensive. This 
means that the functions of variables of state can be 
expand in power series with the small parameter 1 /N 
(N ^> 1) at large finite values of the variables V, z 
for z extensive and large values of V for z intensive. 
Then the extensive function A is said to be the series 
(a > 0) 

A(T,V,z,N) = N[a(T,v,S) + 0(N- a )} 

= Na(T, v, z) z-extensive, (50) 

A(T,V,z,N) = N[a{T,v,z) + 0{N- a )} 

= Na(T, v, z) z-intensive, (51) 

and the intensive function <f> is called the series 

<f>(T,V,z,N) = 0(T,v,z)+O(N- a ) 

= <p(T : v, z) z-extensive, (52) 

c(>(T,V,z,N) = 4>{T,v,z) + 0{N~ a ) 

= 4>(T, v, z) z-intensive, (53) 

where v = V/N is the specific volume, z — z/N is the 
specific z, and a = A/N is the specific A. Comparing 
(|50"|) . (|5T|) and ([55]). (|57|) . it follows that the entropy is 
a homogeneous function of degree 1 only in the ther- 
modynamic limit. Thus the zeroth law of thermo- 
dynamics, the Euler theorem and the Gibbs-Duhem 
relation can be satisfied only in the thermodynamic 
limit. 



B. Ideal gas in canonical ensemble 

The thermodynamic properties of the Renyi statis- 
tics in the canonical ensemble (T, V, z, N) can be in- 
vestigated by the exact analytic functions of the vari- 
ables of state for the classical nonrelativistic ideal gas. 
The main relations for the ideal gas were obtained by 
us in Ref . [2| . Here we express them through the vari- 
able z = l/(q — 1) and study the thermodynamical 
properties of such a system. The Hamiltonian and 
the dimensionless infinitesimal element of phase space 
for the ideal gas of N identical particles can be written 
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The phase space distribution function and the 
norm functions Zr and (H) determined from the sys- 
tem of equations lflU|) and (TlTj) take the form 
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where 
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The value of z in Eq. ([59")) is restricted by the condition 
z < — -AT. The Gibbs partition function in Eq. ([56]) 
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In the following the subscript G indicates the Gibbs 
statistics. In the Gibbs limit z — > ±oo, A" = 
const, V — const the Renyi statistics resembles the 
usual Boltzmann-Gibbs statistics. In this limit the 
phase space distribution function (|55p takes the form 



fcT 2m 



0|z->±oo = Qg = -7T- e 
^G 
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where the norm functions l|56[) and (|57p are 
^flU-±oo = and (iJ> = (H) G . The function (pj, 
d^J in the Gibbs limit is a ZjA r|,^ ±00 = e 3JV/2 . 

The thermodynamic potential of the canonical en- 
semble for the classical ideal gas of the Renyi statistics 
is found to be 
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Here ^p(z) is the psi-function and n G is the chemical 
potential of the Gibbs statistics. It takes the form 
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In the Gibbs limit the chemical potential (|66|) is 
mU-»±oo = MG- The variable X (|2Tj) for the ideal 
gas of the Renyi statistics can be written as 



X = kTX a 
where X a = d(\n.a z ^) / dz and 
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where Fq = —kT In Zq is the free energy of the Gibbs 
statistics. In the Gibbs limit the free energy (|62|) is 
reduced to the form F\ z ^± OQ — Fq- From Eq. (f2Ul) . 
we obtain the entropy S and the pressure p 
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In the Gibbs limit the entropy and pressure can be 
written as S\ z ^±oo — S G and p = p G . The chemical 
potential (f2Tj) for the ideal gas of the Renyi statistics 
in the canonical ensemble can be written as 



u = u G - kT(A^ - -), 



(66) 
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In the Gibbs limit the variable (|T0j) takes the form 
AT|z^±oo = 0. The heat capacity |23|) for the ideal 
gas of the Renyi statistics is now 
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The heat capacity for Renyi statistics is the same as 
for Gibbs statistics. 

The A^-particle distribution function of the classical 
ideal gas in the canonical ensemble for Renyi statistics 
can be defined as 
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which must be normalized to unity 

d 3 p 1 ---d 3 p N f{p 1 ,...,p N ) = 1. (75) 
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In the Gibbs limit the iV-particle distribution function 
([71)1 is reduced to the usual Maxwell-Boltzmann form 
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Integrating Eq. ([74]) with respect to momenta 
p*2 , ■ ■ ■ , pW j we get the reduced single-particle distri- 
bution function 
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In the Gibbs limit we obtain the usual Maxwell- 
Boltzmann one-particle distribution function 
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where 4,ivU-»±oo = 1. 

Let us investigate the zeroth law of thermodynam- 
ics, the principle of additivity and the principle of sta- 
tistical independence on the base of the exact analyt- 
ical results for the ideal gas. Divide the system into 
two parts under the conditions (|2"4"]) - (f2"6"]) . Then from 
Eqs. P> J) -P) ]) and ([65]). ([65]). it follows that the Renyi 
entropy ([63]) and the Gibbs one ([65]) are nonadditive 
in finite systems 
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Comparing Eqs. (EH), ((Ml) and (H2), (J37J) , we see that 
the entropies S and Sq are not homogeneous func- 
tions. Therefore the zeroth law of thermodynamics is 
violated both for the Renyi and Gibbs statistics of the 
finite ideal gas in the canonical ensemble. Moreover, 



the canonical phase space distribution functions for 
the Renyi and Gibbs statistics do not factorize 



8 7^ £>i£>2 



1 



1 



Zr,iZr,2 

1 



(85) 



k,Ti{zi + l 
1 



I.Ni) 



i=i 

N 2 

E 



kT 2 (z 2 + 1 

QG 7^ QG,lQG,2 = 



N 2 ) 



i=l 

2m 



1 



N _ 2 
l "i 

kT 2—1 2m 



Z G .lZ G2 



(86) 



where Zr ^ Zr^Zr^- Thus the principle of statisti- 
cal independence is violated, it is polluted by finite- 
size effects. This completes the proof of the inequali- 
ties ([4"T1) . 

Let us now investigate the Euler theorem for the 
ideal gas in the canonical ensemble. For the Gibbs 
statistics combining Eqs. ([BT]). ([64]) . ([65). and ([69]) . 

we get 

TSg + O g = (H) g +PgV~^ g N, (87) 
O g = kT\}nN\ + N -Nip(N + l)]. 



For the Renyi statistics from Eqs. (|57J). f55|) and 
([72]) . we have in the case of z extensive 



TS + O^ = (H 



O 



pV" + Xz - /xiV, 

= o G + 

+ kT[— lna z jv 



(89) 



+ X a z + VV] (90) 



and in the case of z intensive 



TS + O^ = (H)+pV-tiN, (91) 
G + fcT[- In a a ,jv + VV]. ( 92 ) 



W _ 

R 



Thus the Euler theorem is not valid for the finite 
ideal gas in the canonical ensemble both for the Renyi 
statistics and the Gibbs one. This completes the proof 
of the inequalities 



C. Canonical ideal gas in thermodynamical limit 
with z being intensive 

Consider the ideal gas in the canonical ensemble for 
the Renyi statistics in the thermodynamic limit with 
z intensive: V — > oo, N — > oo and v = V/N = const, 
z = const. For the sake of convenience let us consider 
only the case z > 0. Then the phase space distribution 
function ([55]) and the norm functions ([56"]) , (|6"D|) take 
the form 



1 



1 



N _n 

n 



kT(z 



E 



§iV) frf 2m 



(93) 
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where 



T(z + l) e 



z + l+^N 



V2tt(z 
(gve) N 



1 + 2 

x ~ 2 

mkT 
2nh 2 



~ N 



^G 



and 



a z ,N 



Tjz+T) e z+1 +i N 
^{z + l) z {z + l + lN)h 



(94) 
(95) 

(96) 



In the Gibbs limit z — ► oo and N = const for the dis- 
tribution function (f9"3"|) and the norm function ([Ml) , 
we obtain g| z _,oo = PG and ZrJ^og = Z G , respec- 
tively. Then for Renyi statistics with z intensive the 
Gibbs limit commutes with the thermodynamic limit 
contrary to the Tsallis statistics [22| ■ 

In the thermodynamic limit with z-intensive the 
mean energy of the system ([57)1 and the heat ca- 
pacity ([731) are unchanged, (H) = (3/2)kTN and 
CvzN = (3/2)kN, respectively. However the free en- 
ergy ([62]) is modified 

(97) 
(98) 



F = F G + 0(lnN), 
F G = ~kTN\nZ G , 

O(lnJV) = 2 feTln ( z+ 1 " 



fcTln 



T(z + 1) 



3 -N 

e z+l 



n2n(z + iy 



+ 0(-). (99) 



The entropy (|63[) in the thermodynamical limit with 
z-intensive can be written as 



S 
S G 

OQnN) 



S G + 0(]nN), 



= /,-V \\nZ G + - 



k\n z + 1 + -N 



feln 



r(z + 1) e z+1 



2tt(z + l) 2 



(100) 
(101) 



(102) 



The pressure (I64[) in the thermodynamical limit 
with z-intensive takes the form 



P 



PG 



kT 

v 



(103) 



The chemical potential (166p in the thermodynami- 
cal limit with z-intensive is 



/' 



ji G = —kT(hiZ G — l 



<4> 



kT 



2N 



1 



3 N 
2 z+l 



(104) 
(105) 



0{ W ). (106) 



The variable (|70jl in the thermodynamical limit 
with z-intensive can be written as 



X = X N + O(N°), X o = 0, 



(107) 



O(N ) 



kT[^(z 
kT 



1) 



z + l 



ln(z + 1 
N 



2 V2z + 1 



Oijjs). (108) 



The functions (j97)) - l|107[) were obtained by expand- 
ing them in power series for small 1/N. For exten- 
sive functions the term proportional to N and for in- 
tensive functions the term proportional to N° were 
kept. Note that the chemical potential (| 104|) and the 
thermodynamic variable (|107p are not consistent with 
the thermodynamical relations (f2"Tj) with respect to 
the thermodynamical potential ([9T]) . In the Gibbs 
limit for the thermodynamic variable (|107p . we have 

X\ z ^ oo =0. 

In the thermodynamical limit with z-intensive the 
TV-particle distribution function (fT4"|) of the classical 
ideal gas for Renyi statistics in the canonical ensemble 
is unchanged 



A' 



/(pi, . . . ,p N ) 



h 3N Z R 

1 



(109) 



N rP 

E Pi 
2m 



kT{z + 1 + IN) ^ 



However the single-particle distribution function ([77 
is reduced to the usual Maxwell-Boltzmann form 



f(p) = f G (p) = 



1 



(27rmfcT) 3 / 2 



g 2mkT t 



(110) 



In the Gibbs limit z — > oo and N — const for 
the iV-particle distribution function (| 1 09[) . we have 
/(pi, . . . ,p N )\z^oo = /g(pi, ■ • ■ ,Pn) (ESI). 

Let us also investigate the zeroth law of thermo- 
dynamics, the principle of additivity and the princi- 
ple of statistical independence for the ideal gas in the 
thermodynamical limit with z-intensive. Divide the 
system into two parts under the conditions (|24p. (|26p . 
Then from Eqs. (fTUD|) - (fTUT|) and (S5]), it follows that 
the Renyi entropy (jlOOp and the Gibbs one (jlOip are 
additive (extensive) 



S — Si + S2 = S G ,i 
S G = S G 1 + S G 2, 



s, 



G,2, 



(HI) 

(112) 



where Z G = Z G ^\ = Z G ^- Here terms proportional 
to N were kept only. Therefore for the Renyi and 
Gibbs statistics in the thermodynamical limit with 
z-intensive the zeroth law of thermodynamics is sat- 
isfied. Comparing Eqs. (fTTTj) , (fTT2|) and (36]), (f37|) . 
we see that the entropies S and S G are homogeneous 
functions. The phase space distribution function for 
the Gibbs statistics (|6"Tj) , (|B"6"|) factorizes, however, the 
phase space distribution function for the Renyi statis- 
tics (1851) does not factorize 



Q ^ Q1Q2 = 



1 



1 



(113) 



kT(z + 1 



2—i 9rr 



fiVi) frf 2m 
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fcT(z+l + fiV 2 )^ 



N 2 -2 

V" EL 

2m 



kT 22 ' 



go 



QG,\QG,2 = -z- e 

^G 



(114) 



where 7^ Zr^Zr^ and = Zg,\Zg,2- This 
means that the principle of statistical independence 
is not satisfied for the Renyi statistics at finite N\ 
and N2. 

Combining Eqs. |jS7J|, (ESI), (ED), and ([52]), we obtain 
the Euler theorem 



TSg 



(H)+pV-fj,N, = OQjxN), (115) 

(H) g +PgV-vgN, O g = 0. (116) 



Here only terms proportional to TV are kept. Thus 
for the Renyi and Gibbs statistics in the thermody- 
namical limit with z-intensive the zeroth law of ther- 
modynamics, the principle of additivity and the Euler 
theorem are satisfied. They are valid whenever the 
entropy (|37[) is a homogeneous function in the ther- 
modynamic limit. 



The entropy ([63)) and the pressure (|64[) in the ther- 
modynamical limit with z-extensive are 

S = S G + O(N ) : (122) 

0(iV°) = -* in (l + 1)4.0(1), (123) 

(124) 



V 



Pg 



kT 

v 



The chemical potential ([6"6"]) and the variable ([70)) 
in the thermodynamical limit with z-extensive are 



M = Mg + 0(1), 



(125) 



fcT 



2iV 



' T+0(^) (126) 



and 



X = 



X AT° + O(-), 
fcT 1 



^o = 0, (127) 
O(^). (128) 



D. Ideal gas in thermodynamical limit with z 
being extensive 

Consider now the ideal gas in the canonical ensem- 
ble for the Renyi statistics in the thermodynamic limit 
with z extensive: V — > 00, N 00, z — * ±00 and 
v = V/N = const, z = z/N = const. Then the phase 
space distribution function ()55[) and the norm func- 
tions ([55]). flSSJ), and ([HH]) take the form 



where 



1 



Zj 



a z ,N 



1 



N 



Pi 



-, zN 



kT(z + l)N^2m 



e% N [ 1 



25 



e* w |l + 



25 



zN-1/2 



-1/2 



G ■ 



(117) 

(118) 
(119) 



The partition function Zq is given in Eq. (|95[) . In the 
Gibbs limit 5 — > ±00 and N = const for the distribu- 
tion function (|117p and the norm function (|118p , we 
obtain g\z^±oo = Qg and Z R \z^±< x = Z G , respec- 
tively. Then in Renyi statistics with z extensive the 
Gibbs limit is commutative with the thermodynamic 
limit contrary to the Tsallis statistics [22| . 

The mean energy of the system ([57)) and the heat 
capacity ([73")) are unchanged, (H) = (3/2)kTN and 
CvzN = (3/2)feiV, respectively. The free energy (1521 
in the thermodynamical limit with z extensive can be 
written as 



0(N V ) 



F = F G - 
0- kT 



0(N% 



In 1 



3 

25 



(120) 

+ 0(1). (121) 



Here for the extensive functions (| 1 20[) . I|122[) only 
terms proportional to N are kept, and for the inten- 
sive functions p25[) , (|127p terms proportional to A^ . 
Note that the chemical potential (|125[) obtained in the 
thermodynamic limit does not agree with the ther- 
modynamical relation (|2ip and the thermodynamical 
potential (fT20j) . 

The A^-particle distribution function ([73)) of the 
classical ideal gas for Renyi statistics in the thermo- 
dynamical limit with z extensive takes the form 



f(pi,...,p N ) = 



h 3N Zr- 



1 - 



N 



(129) 



A' 



s2 



kT(z- 



)N ^ 2m 



At the same time the single-particle distribution func- 
tion (JTTJ is reduced to the Maxwcll-Boltzmann form 



m = fern 



(2nmkT) 3 / 2 



(130) 



where d z _N = 1. In the Gibbs limit 5 — ► ±00 and 
N = const for the ./V-particle distribution function 
(|129j) . we have /(pi, . . . ,pn)\z->±oo = /g(pl, • ■ ■ ,Pn) 
(US)). 

Finally, let us investigate the zeroth law of thermo- 
dynamics, the principle of additivity and the princi- 
ple of statistical independence for the ideal gas in the 
thermodynamical limit with z extensive. Divide the 
system into two parts under the conditions (|24|) , (123)) . 
Then from Eqs. ([TTTjl . (fTTH)) and (fl"22)) . (fT23l) . it fol- 
lows that the Renyi entropy (|122p and the Gibbs one 
(|10ip are additive but the distribution function l|117p 
does not factorize 



S — Si + S2 ~ Sg,i + Sq,2 



(131) 
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and 



1 



kT(z + f )N 1 



Ni -a 

^ 2m 



i=i 



kT(z + §)JV 2 



N 2 -2 

2m 



3 = 1 



ZN 2 



(132) 



where the norm function pi8[) factorizes, Zji = 
Zr,iZr,2- The phase space distribution functioni, 
however, does not factorize. Therefore statistical in- 
dependence is not achieved. Combining Eqs. (|122p ~ 
(|127p . we obtain the Euler theorem for z extensive 



TS= (H) + pV + Xz - fiN, 



(133) 



where X = 0. In Eqs. (fT5Tj) , (fH3")l was kept only 
the terms proportional to the first power of N. Thus 
for the Renyi statistics in the thermodynamical limit 
with z extensive the zeroth law of thermodynamics, 
the principle of additivity and the Euler theorem are 
satisfied. They are valid whenever the entropy (|122[) is 
equivalent with (|36[) and it is a homogeneous function 
of the first degree. 

By the exact analytical results for the ideal gas, it 
was proved that the Renyi statistics in the canonical 
ensemble is thermodynamically self-consistent in the 
thermodynamical limit, for z being either extensive or 
intensive. In this case the ideal gas in the canonical en- 
semble for Renyi statistics is equivalent with the ideal 
gas for the Gibbs statistics, except for the N-particle 
phase space distribution function, which on the other 
hand lacks the physical sense of a thermodynamical 
limit. 



III. MICROCANONICAL ENSEMBLE 

(E,V,z,N) 

A. General formalism 

Consider the equilibrium statistical ensemble of the 
closed energetically isolated systems of N particles 
at the constant volume V and the constant thermo- 
dynamic coordinate z. Suppose, these systems have 
identical energy E within AE <C E; then the phase 
space distribution function g be distinct from zero 
only in that region of phase space, D, which is de- 
fined by the inequality E < H < E + AE, with H 
being the Hamiltonian of the system. 

To express the equilibrium phase space distribution 
function through the macroscopic variables of state 
of the microcanonical ensemble we use the thermody- 
namical method based on the fundamental equation 
of thermodynamics [l4| . The fundamental equation of 
thermodynamics at constant values of E, V, z, N can 
be written as 



{dS)E,V,z,N — 0. 



(134) 



Using Eq. {l}, we obtain 

dT—dg+—dz = 0. (135) 

QQ OZ 

D 

Since in the microcanonical ensemble these differen- 
tials vanish, dz — dg — 0, we obtain 



6S_ 

Sg 



dS_ 

Or: 



(136) 



where a\ and ai are certain constants. Substituting 
Eq. ((T|) for the Renyi entropy in the first equation of 
(|136[) and in its expectation value, we get 

= -k(z+l)Q 1 ' z e& =ai, (137) 

og 



/$ S 
dTg— = — k(z+l) = ai. 
dg 



(138) 



Combining Eqs. (fHT)) . [fTggj) . and ©, we obtain the 
phase space distribution function and the statistical 
weight for the Renyi statistics in the framework of 
the classical statistical mechanics [2[ 

1 



where 



W 



J dT = J dTA(H-E). 

D 



(139) 



(140) 



Then the Renyi entropy is reduced to the familiar 
expression 

S = k\nW = S Gl (141) 

where Sq is the Gibbs entropy. Since the Renyi en- 
tropy (I14ip is a thermodynamical potential of the mi- 
crocanonical ensemble and it is equivalent with the 
Gibbs entropy, and it is independent of the variable 
of the state z, we obtain that in the microcanonical 
ensemble the Renyi statistics is equivalent with the 
usual Boltzmann-Gibbs statistics and it is indepen- 
dent of the variable of the state z (a2 = 0). Therefore 
in the microcanonical ensemble we have only three 
variables of state (E, V, N). From Eq. 1(159]) it follows 
that the expectation value of a dynamical variable A 
([3J can be written as 



{A) = J dTgA = ±J dTA(H - E)A, 

D 



(142) 



where (A) = (A)g- For the distribution function (|139|) 
and the entropy (|141[) . we have the fundamental equa- 
tion of thermodynamics 



where 



TdS = dE+ pdV - ndN, 
I - f^.\ - _L 

f - [dE) VN -%' 



p = T 



(i = -T 



dV 



PG, 



E.N 



dS 
dN 



(143) 

(144) 
(145) 
(146) 



E.V 
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Here (H) = E and X = T(dS/dz) = 0. 

The fundamental equation of thermodynamics 
(|143p provides the first and second laws of thermo- 
dynamics 

6Q = TdS, SQ = d(H) + pdV - fjdN. (147) 

The heat capacity, SQ — CdT, in the microcanonical 
ensemble at the constant values of E, V, z, N can be 
written as 



C 



1 



EVN 



(dT\ 
\dEJ 



V,N 



T 2(d 2 S\ 

± \dW)v,N 



(148) 



where Cevn 

Let us prove the zeroth law of thermodynamics for 
the Renyi statistics in the microcanonical ensemble. 
Divide the system into two parts, H = Hi + H 2 , under 
the conditions 

E = Et =E 2 , N = N 1 + N 2 , V = V 1 +V 2 . (149) 

The heat transfer SQ is additive. Therefore, from the 
second law of thermodynamics p47[) it follows that 
the zeroth law, T = T\ = T 2 , is valid whenever the 
total differential of entropy is additive dS = dSi + dS 2 
(see Eqs. (J27J), (|25j)l. Let us show that Eq. (J2SJ) is valid 
only if the entropy S is additive function 

S(E,V,N) = S 1 (E 1 ,V 1 ,N 1 ) + S 2 (E 2 ,V 2 ,N 2 ). (150) 



Differentiating Eq. (|150p with respect to variables of 
state E, V and N, we obtain 



as 


dSi 


8S 2 




OE 


dEi 


dE 2 ' 




dS 


dSi 


dS 2 




dV 


dVi 


~ dV 2 ' 




dS 


dSi 


dS 2 




dN 


dNi 


~ 0N 2 ' 




T = 


Ti = 


--T 2 , 


E 


P = 


Pi = 


-P2, 


V 


P = 


Mi = 


-- M2, 


N 



E = E X - 
V = V 1 + 
N = Ni 



-E 2 
V 2 , 



-E 2 , 
V 2 , 
f N 2 . 



(151) 
(152) 
N 2 (153) 



(154) 
(155) 
(156) 



Then from Eqs. lfi~ST |) -([l~53" j) and the differentials of 
Eq. <JHHD, we get 



dS = dS!+dS 2 . 



(157) 



Therefore the zeroth law of thermodynamics, T = 
Ti = T 2 , in the microcanonical ensemble for the Renyi 
statistics is satisfied whenever the entropy is additive, 
S — S\ + S 2 . 

Let us show the correspondence in the microcanon- 
ical ensemble between the additivity of entropy (|150[) 
and its homogeneity. Consider the Renyi (Gibbs) en- 
tropy in the microcanonical ensemble be a first order 
homogeneous function of variables E, V, N 

S(E,XV,XN) = \S(E,V,N), A = — . (158) 



Divide the system into two parts under the conditions 
For A, we have Eq. fgg]). Then the intensive 
specific variables correspond to the additive (exten- 
sive) variables 



XE — X\E\ — X 2 E 2 



XV = XiVi = X 2 V 2l (159) 



where Ai = 1/iVi and A2 = 1/N 2 . Under the condi- 
tions of Eqs. (fTSTij) , (TT55|) . and (TT591) . ^ follows that 
the specific entropy AS* is intensive function 

XS{E,V,N) = XiSi(E 1 ,Vi,N 1 ) = 

= X 2 S 2 (E 2 ,V 2 ,N 2 ). (160) 



Conversely, from Eqs. (jlSO]) . ([155]) . and |TSS]), we ob- 
tain the additivity of entropy (|150p . Thus the zeroth 
law of thermodynamics for Renyi statistics in the mi- 
crocanonical ensemble is valid whenever the entropy 
(or the thermodynamic potential) is a homogeneous 
function of first degree of the extensive variables of 
state Him. 

If the entropy is a homogeneous function of de- 
gree 1 (|158[) , then S satisfies the Euler theorem and 
the Gibbs-Duhem relation. Differentiating Eqs. (| 1 58[) 
with respect to A, and putting A = 1, we obtain the 
well-known Euler theorem for the homogeneous func- 
tions: 



El'—} v(— 

'VdE) VN + \dV /ES 



N 



as 

dN 



E,V 



s. 

(161) 



Using the thermodynamic relations (|144|) - (|146|) for 
the isolated thermodynamic system (E, V, N), we get 
the Euler theorem [29j |: 



TS = E +pV - fiN. 



(162) 



Applying the differential operator with respect to the 
variables of state (E, V, N) on Eq. (|162p . we obtain 
the fundamental equation of thermodynamics (I143[) 
and the Gibbs-Duhem relation [29j 



SdT = Vdp - Ndfi. 



(163) 



Equation (|163p means that the variables T, fi and p 
are not independent. The fundamental equation of 
thermodynamics (|143p provides the first and the sec- 
ond laws of thermodynamics (|147p . Therefore in the 
microcanonical ensemble the Euler theorem (|162p and 
the Gibbs-Duhem relation p63p result from the homo- 
geneity property of the entropy (|158p . 

Dividing the finite system of the microcanonical en- 
semble into two dynamically independent subsystems, 
H = Hi + H 2 , under the conditions (| 149)1 , we have a 
convolution of the statistical weight W 



W(E, V, N) 



Ni\N 2 \ V^Y?* 
N\ V N 



XL 

x J dEiW^EuVuN^W^E - E U V 2 ,N 2 ), (164) 



where Wi = J A(Hi — Ei)dTi. Hence the statisti- 
cal weight (|140p and the distribution function (|139[) 
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do not factorize, W ^ W\Wi and g ^ Q1Q2, respec- 
tively. Therefore the Renyi (Gibbs) entropy (| 141 [) is 
nonadditive (nonextensive) function and the proper 
temperature (|144p is not an intensive variable 



S^S 1 + S 2 , T^T^T 2 



(165) 



Thus the Gibbs and Renyi statistics for finite sys- 
tems violate the zeroth law of thermodynamics (see 
Eq. (|165[) ). Moreover, the Euler theorem is also not 
valid 



TS y£ (H) + P V- fiN. 



(166) 



For finite values of the variables of state E, V, N the 
thcrmodynamical potential (|141[) of the microcanon- 
ical ensemble is not a homogenous function of first 
order. The proof of inequalities p65[) . (|166p is given 
below on the base of the ideal gas. 

In the general case, the Renyi (Gibbs) entropy (|14ip 
of the variables of state is not a homogeneous func- 
tion. It can become homogeneous only in the ther- 
modynamic limit, iV — > 00, E — > oc, V — > 00 and 
v = V/N = const, z = z/N — const. This means that 
the functions of variables of state can be expand in 
power series with the small parameter 1/N (N ^> 1) 
at large finite values of the variables E, V. Then the 
extensive variables A and the intensive ones <p can be 
written (a > 0) as 

A(E,V,N) = N[a{e,v) + 0{N' a )} 

= Na(e,v), (167) 

4>{E,V,N) = <f>(e,v) + 0(N- a ) 

<f>(e,v). (168) 

where e = E/N is the specific energy, v = V/N is the 
specific volume, and a = A/N is the specific A. Com- 
paring (|167p and (|158p , it follows that the entropy is a 
homogeneous function of degree 1 only in the thermo- 
dynamic limit. Thus the zeroth law of thermodynam- 
ics, the Euler theorem and the Gibbs-Duhem relation 
for the Renyi (Gibbs) statistics can be satisfied only 
in the thermodynamic limit. 



B. Ideal gas in microcanonical ensemble 

The thermodynamic properties of the Renyi 
(Gibbs) statistics in the microcanonical ensemble 
(E, V, N) can be investigated by the exact solved an- 
alytical functions of the variables of state for the clas- 
sical nonrelativistic ideal gas. The main relations for 
the ideal gas of the Gibbs statistics can be found 
in [l4j]. The phase space distribution function (|139p 
and the statistical weight (| 140[) are given by [30J 



Q - 



1 




d 3 Pl . . . d 3 p N 



V N ( m 



\" E lN-l 

TON)' 



(169) 



where m is the particle's rest mass. The entropy (|141l) 
and the temperature (|144[) take the form 



S = Sq = k In 



T = T r = 



yN 

E 



r(fAT) 



UN -l)k 



, (170) 
(171) 



The pressure (11451) . the chemical potential Q145P and 
the heat capacity (|148|) can be written as 



PG 
kT 



*kT 
V ' 



-fcTln 



V 



mE 
2nh 2 



3/2' 



V(JV + i) + |v 



c 



EVN 



C G 



EVN 



-N — I] k 



E 
T' 



(172) 

(173) 
(174) 



Divide the system into two parts under the condi- 
tions (HI]) . Then from Eqs. (T6l? H l"71" 1) . it follows 
that the Renyi (Gibbs) entropy (|170p is nonadditive, 
the temperature (|17ip is nonintensive and the phase 
space distribution function does not factorize 

S £ Si+Sn = k\n.(WiW 2 ), T^T X #T 2 ,(175) 

g * 9192 = w^w 2 > (176) 

where W ^ W1W2. Using Eqs. (fT70| - (fT73| . we have 

(177) 

N \ 

Or = Ell 



TS + R = E+pV-fiN, 
N 



In -i 



EN 



2 iv 



1 



iP(N+l) + ^N) 



.(178) 



Thus the zeroth law of thermodynamics, the principle 
of additivity and the Euler theorem are not satisfied 
in the microcanonical ensemble for the finite ideal gas 
of the Renyi (Gibbs) statistics. This completes the 
proof of the inequalities (|165p , (|166p . 

Consider now the ideal gas in the microcanoni- 
cal ensemble for the Renyi statistics in the thermo- 
dynamic limit E ^ oo, V — > 00, N 00 and 
e = E/N = const, v = V/N = const. Then the 
phase space distribution function and the statistical 
weight p69p can be written as 



Q= W , W=,r 



N 



mee 



3/2 



(179) 



The entropy (|14ip . the temperature (|144p . the pres- 
sure (|145p . the chemical potential (|145[) and the heat 
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capacity (|148[) in the thermodynamic limit are found 
to be 



S = Sq = kNhiw, 
3 ft' 

kT 2e 



T = T G = 



v 3 v ' 



P = Pg — 



fx = hg = —kT lnw 



C 



EVN 



= C<i VN = ~kN = N^. 



2 J ' 



(180) 
(181) 

(182) 

(183) 

(184) 



In the thermodynamical limit the Euler theorem 
(fTTTD . (fTTSD is valid 



(185) 



where s = S/N . 

It is easy to prove that in the thermodynamic limit 
under the conditions E = E\ + E2, V = V\ + V% and 
N = Ni + N-2 or e = E\ = e 2 , v = v\ = v 2 the Renyi 
entropy (|180p is additive, the phase space distribution 
function (|179[) factorizes and the temperature (| 18 1 [) is 
intensive, i.e., the zeroth law of thermodynamics is 
valid 



S = Si + S 2 



T = T 1= T 2 



(186) 



Here g = qiq 2 and W = W\W 2 . Therefore the 
Renyi (Gibbs) statistics in the microcanonical ensem- 
ble completely satisfies all requirements of the equilib- 
rium thermodynamics. Let us note that the zeroth law 
of thermodynamics, the principle of additivity and the 
Euler theorem for the Renyi (Gibbs) statistics in the 
microcanonical ensemble are satisfied due to the fact 
that in the thermodynamic limit the thermodynamical 
potential of the microcanonical ensemble, the entropy 
(|180|) . becomes a homogeneous function of first degree 
of the extensive variables of state. 



IV. CONCLUSIONS 

To conclude, this paper has explored the Renyi 
statistics in the canonical and microcanonical ensem- 
bles in the general case and for the classical ideal gas. 
The exact analytical results for the ideal gas were ob- 
tained both for a finite system and in the thermody- 
namic limit. In the canonical ensemble the ideal gas 
in a thermodynamic limit was studied in two cases: 
when the variable z is extensive and when it is in- 
tensive. The problem of the connection between the 
statistical mechanics based on the Renyi entropy and 
the equilibrium thermodynamics in the canonical and 
microcanonical ensembles was analyzed. 

Here we summarize the main results of this study. 
The phase space distribution functions for the canoni- 
cal and microcanonical ensembles were derived by the 
method based on the fundamental equation of thermo- 
dynamics instead of the Jaynes principle. Both these 



methods give the same results. For the Renyi statis- 
tics in the canonical ensemble we obtain the power-law 
phase space distribution function fixed by two norm 
functions instead of one partition function given in 
the Gibbs statistics. However in the microcanonical 
ensemble we have the usual equiprobability distribu- 
tion function. The main thermodynamical relations 
and laws for a statistical ensemble of the Renyi statis- 
tics were obtained from the ensemble averages and 
the phase space distribution function. The connec- 
tion of the Renyi statistics and the equilibrium ther- 
modynamics was established. By the exact analytical 
results for the ideal gas, it was shown that the zeroth 
law of thermodynamics, the principle of additivity, the 
Euler theorem and the Gibbs-Duhem relation are valid 
only if the thermodynamical potential of the statisti- 
cal ensemble is a homogeneous function of degree 1 
of its extensive variables of state. For finite systems 
the thermodynamical potential is an inhomogeneous 
function and the Renyi (Gibbs) statistics is neither 
in the canonical nor in the microcanonical ensemble 
is thermodynamically self-consistent. The thermody- 
namical potential becomes a homogenous function of 
first degree only in the thermodynamic limit. In this 
case it is additive, and the Renyi (Gibbs) statistics is 
thermodynamically self-consistent. 

It was further shown that the thermodynamical po- 
tential of the microcanonical ensemble, the Renyi en- 
tropy, is equivalent with the Gibbs entropy and is 
independent of the variable of state z. Therefore 
in the microcanonical ensemble the Renyi statistics 
completely coincides with the usual Boltzmann-Gibbs 
statistics. It was found that for finite systems in the 
canonical ensemble the entropic index z must be an 
extensive variable of state for the Euler theorem be- 
ing consistent with the fundamental equation of ther- 
modynamics and the Gibbs-Duhem relation. It was 
revealed that in the canonical ensemble in the ther- 
modynamical limit both for z extensive and z inten- 
sive the ideal gas of the Renyi statistics is equiva- 
lent with the ideal gas of the Gibbs statistics, except 
for the phase space distribution function, which nev- 
ertheless in the thermodynamic limit lacks physical 
sense. Evidently, the Renyi statistics in the canon- 
ical and microcanonical ensembles is thermodynam- 
ically self-consistent whenever it coincides with the 
usual Boltzmann-Gibbs statistics. Therefore we con- 
clude that the Renyi statistics duplicates the ther- 
modynamic relations of the familiar Boltzmann-Gibbs 
statistics in the thermodynamic limit, and leads to no 
new result in this respect. This conclusion relies on 
the particular approach, that the extra parameter in 
the entropy formula, q is treated as related to a vari- 
able of state, z = l/(q — 1). 
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